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ABSTRACT 

We present the size evolution of passively evolving galaxies at z ~ 2 identified in Wide-Field Camera 3 imaging 
from the Early Release Science program. Our sample was constructed using an analog to the passive BzK galaxy 
selection criterion, which isolates galaxies with little or no ongoing star formation at z > 1.5. We identify 30 
galaxies in ~40 arcmin 2 to // < 25 mag. By fitting the 10-band Hubble Space Telescope photometry from 
0.22 /x m < l 0 bs <1.6 /x m with stellar population synthesis models, we simultaneously determine photometric 
redshift, stellar mass, and a bevy of other population parameters. Based on the six galaxies with published 
spectroscopic redshifts, we estimate a typical redshift uncertainty of ^0.033( 1 +z). We determine effective radii from 
Sersic profile fits to the //-band image using an empirical point-spread function. By supplementing our data with 
published samples, we propose a mass-dependent size evolution model for passively evolving galaxies, where the 
most massive galaxies (M* ~ 10 1 1 M G ) undergo the strongest evolution from z ~ 2 to the present. Parameterizing 
the size evolution as (1 + z) _ “, we find a tentative scaling of a ss (—0.6 ± 0.7) + (0.9 ± 0.4) log(M*/10 9 M 0 ), 
where the relatively large uncertainties reflect the poor sampling in stellar mass due to the low numbers of high- 
redshift systems. We discuss the implications of this result for the redshift evolution of the M*-R e relation for red 
galaxies. 

Key words: galaxies: evolution - galaxies: fundamental parameters - galaxies: structure 
Online-only material: color figures 


1. INTRODUCTION 

Understanding the red, passively evolving galaxies at in- 
termediate redshifts (z ~ 2) is one of the outstanding chal- 
lenges of galaxy evolution studies. Early expectations of a high- 
luminosity phase associated with the formation of spheroids 
on a free-fall timescale (Eggen et al. 1962) have long given 
way to a framework in which the spheroids are assembled over 
an extended period of time. The identification of fairly high 


* Based on observations made with the NASA/ESA Hubble Space Telescope, 
obtained from the Data Archive at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 


space densities of massive and passively evolving galaxies at 
Z > 1.5 (e.g., Glazebrook et al. 2004; Cimatti et al. 2004) re- 
vealed weaknesses in the early semi-analytic models of their 
formation. Large-scale surveys, both locally and out to z ~ 1, 
provide a fairly clear view of the evolution in the number density 
of red sequence galaxies and the global stellar mass density in 
passive systems (e.g.. Bell et al. 2006; Faber et al. 2007; Brown 
et al. 2003, 2007). 

Deep spectroscopic studies of red galaxies (e.g., McCarthy 
et al. 2004; Cimatti et al. 2006, 2008) estimated stellar ages 
consistent with early formation redshifts (Zf 0 rm < 4), potentially 
at odds with the rapid evolution in stellar mass density at these 
epochs (e.g., Rudnick et al. 2003). It is now clear that the red 


1 


The Astrophysical Journal, 749:53 U lpp), 2012 April 10 


Ryan et al. 


sequence was in place at fairly early epochs (Bell et al. 2004; 
Bremer et al. 2006; Kriek et al. 2008), but evolved strongly in 
the 1 < z < 3 era (e.g., Demarco et al. 2010) and likely since 
Z ~ 1 as well (e.g., Faber et al. 2007). The emerging picture in 
which massive galaxies are assembled via mergers followed by a 
rapid quenching of star formation addresses many of the salient 
properties of the red sequence galaxies (Faber et al. 2007). 

In light of these results it was quite surprising that high 
spatial resolution studies showed that the passive red galaxies 
at intermediate redshifts are systematically smaller than their 
likely present-day counterparts (e.g., Daddi et al. 2005; Trujillo 
et al. 2006; Toft et al. 2007). Observations of 1 < z < 2 passive 
galaxies with either the Advanced Camera of Surveys (ACS) or 
the Near Infrared Camera and Multi-Object Spectrometer on the 
Hubble Space Telescope ( HST) implied sizes that are a factor of 
~2 smaller than equal mass red galaxies today (e.g., Longhetti 
et al. 2007; Cimatti et al. 2008; Damjanov et al. 2009), while 
observations of red galaxies at z > 2 suggest even more dramatic 
evolution (e.g., van Dokkum et al. 2008; Saracco et al. 2010). 

The primary interest in determining the characteristic sizes of 
passive galaxies is the insight that it provides into the dynamical 
state of the system. As one projection of the fundamental plane 
(Djorgovski & Davis 1987; Dressier et al. 1987), the Kormendy 
relation (Kormendy 1977) provides a probe of the dynamics 
of hot stellar systems. The Kormendy diagram for the z > 1 
red sequence galaxies is nearly as tight as it is locally, but is 
offset in both surface brightness and size (e.g., Damjanov et al. 
2009), even when the effects of passive fading of the stellar 
populations are taken into account. The compact sizes and high 
surface brightness of the red galaxies imply stellar densities 
in the centers of these galaxies that are two to three orders 
of magnitude higher than present-day massive red galaxies (e.g., 
van Dokkum et al. 2008; Damjanov et al. 2009). 

The challenge in understanding these results stems from the 
apparent conflict between the requirement for a factor of ~3 
growth in size during an epoch in which the stellar masses, 
population ages, and overall morphologies show little or no 
evolution. A number of models are under active discussion 
(e.g., Naab et al. 2006; Fan et al. 2008; Damjanov et al. 

2009) , but most involve large numbers of late-stage minor 
mergers that grow the galaxies in size without contributing 
large amounts of mass. While it is difficult to conclusively 
rule out various measurement biases which might give rise 
to this apparent trend (e.g., underestimated effective radii or 
overestimated masses), most interpretations are astrophysical in 
origin: (major or minor) merging occurring at different phases 
(e.g., Naab et al. 2006), or expansion due to a significant 
mass loss, either by an active galactic nucleus (Fan et al. 
2008) or stellar winds (Damjanov et al. 2009). Recently, 
Flopkins et al. (2010) presented a semi-analytic model based 
on high-resolution hydrodynamic simulations (Cox et al. 2006), 
which incorporates several astrophysical and observational 
mechanisms for this trend, and concluded that for galaxies with 
a stellar mass of M* Js 1 0 1 1 M G , the late-stage minor merging is 
the dominant mechanism, accounting for ~50% of the apparent 
size evolution. Certainly many of these high-redshift systems do 
in fact show merger-like features (e.g., van Dokkum & Brammer 

2010) or have (multiple) companions (e.g.. Bell et al. 2006), 
and so merging is an attractive explanation for the observed size 
evolution. 

The merging scenario implies that some non-negligible frac- 
tion (<10%; Hopkins et al. 2009) of compact galaxies should 
remain at z ~ 0. However, Trujillo et al. (2009) find that only 


0.03% (= 48/152, 083) of galaxies at z~0.2 have stellar densi- 
ties comparable to these high-redshift galaxies, arguing against 
the merger scenario. To further confound the issue, these galax- 
ies are generally young (~3 Gyr old), giving an approximate 
formation redshift of Zf 0rm ~ 0.3, suggesting that they are not the 
relics of the early universe, but formed from gas-rich, recently 
merged disks. Conversely, Saracco et al. (2010) find that ~62% 
(= 21/34) of their galaxies at 1 < z spe c < 2 are within ler of 
the local R e -M * relation (Shen et al. 2003). Based on this find- 
ing, they suggest that the compact high-redshift galaxies are not 
the progenitors of large field galaxies, but of compact brightest 
cluster galaxies. The absence of compact early-type galaxies at 
z~0 (e.g., Taylor et al. 2010) seems to suggest that the simple 
picture of galaxy merging with some mass-loss expansion is 
not completely correct, in stark contrast to the simulated results 
(e.g., Naab et al. 2009; Hopkins et al. 2010). In this paper, we 
present results from the Early Release Science (ERS) program 
with Wide-Field Camera 3 (WFC3) which suggest that the size 
evolution from z ~ 2 to the present depends on stellar mass. This 
finding should provide an additional observational constraint to 
further refine the theoretical models. 

This paper is organized as follows. In Section 2, we describe 
the observations and ancillary data, in Section 3 we define our 
sample, in Section 4 we detail our stellar population modeling 
and morphological analysis, in Section 5 we describe our mass- 
dependent size evolution model, in Section 6 we discuss several 
key results, and in Section 7 we give a brief summary of this 
work with thoughts for future surveys. Throughout this paper, 
we assume a ACDM cosmology with Qo = 0.27, Qa = 0.73, 
and Hq = 73 km s~ 1 Mpc” 1 (Spergel et al. 2003) and will quote 
all magnitudes in the AB„ magnitude system (Oke & Gunn 
1983). 

2. OBSERVATIONS 

2.1. The Wide-Field Camera 3 Early Release Science 

For this work, we analyze the ERS (PropID: 11359, PI: R. 
W. O’Connell) observations conducted with the WFC3 recently 
installed on the HST. This field covers roughly the northern 
40 arcmin 2 of the Great Observatories Origins Deep Survey 
southern field (GOODS-S; Giavalisco et al. 2004), and leverages 
the existing ACS optical data with near-ultraviolet (NUV) and 
near-infrared (NIR) data with equivalent space-based imaging. 
The main details regarding the data collection, reduction, 
calibration, and mosaicking are presented by Windhorst et al. 
(2011). We briefly summarize the observational aspects critical 
to this work. 

The ERS program utilizes the two complementary modes of 
the WFC3: UVIS and IR channels. The data were collected 
between 2009 September and November and form a 2 x 4 
and 2x5 mosaic, respectively. The final data set consists of 
~40 arcmin 2 of 10-band HST imaging covering 2250 A to 
1.6 pm in wavelength. The v2.0 GOODS high-level science 
products had an original pixel scale of 30 mas and were 3x3 
block summed to produce science and weight maps of equal 
pixel scale to the WFC3 observations (see Windhorst et al. 2011, 
for a more detailed discussion on the rebinning process and its 
motivation). Therefore, our final HST mosaics have 90 mas 
pixels. We refer hereafter to the 10-band WFC3 and ACS filter 
set in the F225W, F275W, F336W, F435W, F606W, F775W, 
F850LP, F098M, F125W, and F160W as UiU 2 U 3 BVi'z'Y s J H 
to simplify the bandpass notation. 
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2.2. Source Catalogs 

We created source catalogs using SExtractor (Bertin & 
Arnouts 1 996) with the H- band mosaic as the detection image 
and each of the HST mosaics as measurement images. We use 
the weight maps produced by MultiDrizzle (Koekemoer et al. 
2002), modified to account for correlated pixel noise (Dickinson 
et al. 2004). For object detection, we require a minimum of 
5 connected pixels each greater than 0.6a above the local 
background and apply a 3 pixel Gaussian smoothing filter. 
For deblending, we adopt a contrast parameter of 1 0 4 with 
64 sub-thresholds. For photometry, we adopt the MAG_AUT0 
measurements with a Kron factor of 2.5, and minimum object 
radius of 3.5 pixels, which reliably recover total fluxes to within 
~6% (Bertin & Arnouts 1996). We use the AB zero points 
given by Kalirai et al. (2009a, 2009b) for the UV1S and IR data, 
respectively. 

3. SAMPLE SELECTION 

We select our passively evolving galaxies using a variation 
on the standard p BzK selection (Daddi et al. 2004), designed to 
identify galaxies at z ~ 2 with little ongoing star formation and 
stellar population ages of ~1 Gyr (Daddi et al. 2005). However, 
to take full advantage of our high spatial resolution and signif- 
icantly deeper ERS data set, and avoid issues associated with 
cataloging images with markedly different angular resolution, 
we use our H - band imaging instead of the ground-based A'-band 
images. Therefore, we modify the classical pBzK selection cri- 
teria to an equivalent p BzH scheme: 

(z — H) — {B — z') < -0.2 - ((// - K)) mag, (1) 

(z — H) > 2.5 - ((// - K )) mag. (2) 

As these p BzK galaxies have maximally old, passively evolv- 
ing systems, we derive a typical color of ((// — K)) = 0.7 mag 
based on a stellar population with an instantaneous star for- 
mation history formed at Zf orm = 10 based on S. Chariot & 
G. Bruzual (2007, in preparation) and Bruzual (2007, hereafter 
CB07). In Figure 1 , we show the (B—z') and (. z'—H ) colors with 
these criteria illustrated as a shaded polygon. To rule out any 
potential image artifacts associated with the ERS field edges, 
we require the objects to be in the portion of the //-band mosaic 
which received the full complement of two orbits per pointing. 
We restrict the H - band magnitude to a relatively conservative 
limit of H ^ 25 mag to ensure a complete and reliable sample 
(Windhorst et al. 2011). With these requirements, we identify 
a final sample of 30 galaxies in 43.1 arcmin 2 and present their 
photometry (Table 1). Of these galaxies, 15 objects have very 
little (or no) flux in the B band. To highlight their red colors and 
visual morphologies, we show a 4"5, three-color stamp of each 
object in Figure 2. 

We show the source counts of BzH galaxies in Figure 3. 
The counts plateau around H — 22.5 mag, which is roughly 
4 mag brighter than the //-band completeness limit (Windhorst 
et al. 201 1). Furthermore, since the GOODS z'-band imaging is 
complete to z! ~ 27 . 3 mag (Dahlen et al. 20 1 0) and our effective 
color limit is (z! — H) ^ 1.8 mag, this plateau is unlikely due 
to incompleteness in the z'-band data, suggesting that the faint 
end of their luminosity function is relatively flat or declining 
at low luminosities. Stutz et al. (2008) note a similar dearth of 
low-mass systems (<2x 10 l() M 0 ), further supporting a flat or 
declining faint-end slope. At the bright end (H < 22 mag), our 



Figure 1. BzH object selection. The small dots represent all objects detected 
and measured by SExtractor. The gray area is the color selection region 
defined by Equations (1) and (2) with a typical color of (( H — K )) = 0.7 mag. 
We identify 30 galaxies strictly meeting the BzH color criteria, for objects 
which were undetected in the B band, we show the 5<r point-source upper limits 
on the ( B — z ! ) color (there are three objects off the plot). The few objects in 
the shaded selection region which are not in our sample either fall below our 
flux limit of H < 25 mag or did not receive the full exposure time. We show 
a maximally old population with an instantaneous burst of star formation and 
the canonical (Coleman et al. 1980) elliptical galaxy template as solid red and 
dashed blue lines, respectively. We indicate redshift tick marks and values in 
the corresponding color. The arrow in the lower right comer shows the expected 
reddening for Ay = 0.5 mag in the rest frame at z = 1-6 for the Calzetti 
et al. (1994) dust law. We show the positions of Galactic stars from the Pickles 
(1998) library as open triangles and the color represents the luminosity class (I: 
magenta, II: red. III: green, IV: blue, V: cyan). 

(A color version of this figure is available in the online journal.) 

survey and modified p BzK selection method produces source 
counts consistent with Lane et al. (2007), based on the UKIRT 
Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS). 
While the UKIDSS/UDS survey is much wider (~0.6 deg 2 ), the 
ERS data push > 2 AB mag deeper into a regime not routinely 
possible from the ground. 

4. ANALYSIS 

4.1. Surface Brightness Models 

To determine the rest-frame optical morphologies of our 
galaxies, we model the two-dimensional light distribution in 
the H band using GalFit (Peng et al. 2002). We determine an 
empirical point-spread function (PSF) from a median stack of 
75 stars identified in the field, based on their full width at half- 
maximum (FWHM) and brightness determined by SExtractor 
(Windhorst et al. 201 1). We fit the standard Sersic profile: 

X(r) = I. e e- b ^ r/r ‘ )ll "~ I \ (3) 

where r e is the effective radius, 23 X,. is the surface brightness 
at the effective radius, n is the Sersic index, and b„ is a 
constant found by numerically solving T(2 n) = 2y(2n. b„) 
with T(x) and y(x, y ) representing the gamma and incomplete 
gamma functions, respectively. Since GalFit is minimizing a 
goodness-of-fit statistic, modeling the uncertainty in each pixel 
is critical for meaningful estimate of the uncertainty in each of 
the fit parameters. Therefore, we transform the scaled weight 
maps produced by MultiDrizzle into uncertainty maps and 

23 We refer to the effective radius as r e when in angular units and as R e when 
converted to physical units. 
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ID 


R.A. a 

/h m 

Decl. a 

( ° , n ) 

B 

(mag) 

V 

(mag) 

V 

(mag) 

t! 

(mag) 

Y s 

(mag) 

J 

(mag) 

H 

(mag) 

Notes 

408 

3 

32 28.02 

-27 40 31.2 

27.52 ± 0.49 

26.18 ±0.12 

25.18 ±0.08 

24.40 ± 0.05 

23.77 ± 0.03 

22.87 ±0.01 

22.54 ±0.01 


606 

3 

32 23.27 

-27 40 45.8 

>28.11 

27.60 ± 0.53 

26.12 ±0.23 

25.36 ±0.14 

25.16 ±0.17 

23.66 ±0.03 

23.10 ±0.02 


1696 

3 

32 42.08 

-27 41 41.2 

>27.36 

>27.29 

26.69 ± 1.12 

26.16 ±0.77 

25.47 ± 0.32 

24.98 ±0.14 

24.08 ± 0.08 


2227 

3 

32 42.34 

-27 42 04.0 

>27.55 

26.85 ± 0.63 

24.53 ±0.13 

24.50 ±0.14 

23.86 ±0.07 

22.91 ±0.02 

22.32 ±0.01 


2377 

3 

32 25.03 

-27 42 09.6 

>28.10 

27.80 ± 0.67 

27.68 ± 1.04 

27.43 ± 1.00 

26.44 ± 0.42 

25.42 ±0.11 

24.52 ± 0.06 

X-ray ID 226 b 

2749 

3 

32 14.92 

-27 42 21.9 

>27.23 

26.40 ± 0.44 

24.92 ±0.19 

24.64 ±0.19 

24.67 ±0.19 

23.22 ± 0.03 

22.69 ± 0.03 


2750 

3 

32 14.88 

-27 42 23.3 

>27.26 

25.70 ±0.22 

24.56 ±0.14 

23.92 ± 0.09 

23.95 ± 0.09 

22.46 ± 0.02 

21.97 ±0.01 

X-ray ID 145 b 

2871 

3 

32 31.09 

-27 42 26.6 

>28.23 

26.85 ± 0.22 

25.78 ±0.14 

24.91 ± 0.08 

24.17 ±0.05 

23.39 ±0.02 

23.07 ±0.02 


3000 

3 

32 43.93 

-27 42 32.4 

26.23 ± 0.47 

25.12 ±0.13 

23.61 ±0.06 

22.64 ± 0.03 

21.99 ±0.02 

21.07 ±0.00 

20.59 ± 0.00 


3152 

3 

32 26.05 

-27 42 36.6 

>27.96 

27.21 ± 0.44 

24.93 ± 0.09 

23.93 ± 0.05 

23.31 ±0.02 

22.15 ±0.01 

21.66 ±0.00 


3237 

3 

32 35.91 

-27 42 40.9 

>28.04 

26.77 ± 0.30 

25.24 ±0.13 

24.60 ± 0.08 

23.90 ±0.05 

22.67 ±0.01 

22.06 ± 0.01 


3360 

3 

32 30.58 

-27 42 43.4 

>28.06 

27.27 ± 0.44 

25.90 ±0.22 

24.65 ± 0.08 

24. 1 1 ± 0.05 

22.78 ±0.01 

22.22 ±0.01 


3376 

3 

32 35.13 

-27 42 37.0 

>28.04 

27.95 ± 0.74 

26.80 ±0.45 

26.69 ± 0.48 

25.40 ±0.16 

24.91 ± 0.07 

24.55 ± 0.06 


3471 

3 

32 27.94 

-27 42 45.7 

28.29 ± 1.34 

26.29 ±0.17 

24.75 ± 0.07 

24.02 ± 0.04 

23.48 ± 0.02 

22.35 ±0.01 

21.91 ±0.00 


3488 

3 

32 25.65 

-27 42 46.8 

>27.97 

>28.20 

>27.63 

26.95 ± 0.70 

26.01 ±0.31 

25.38 ±0.12 

24.93 ±0.10 


3551 

3 

32 36.28 

-27 42 49.4 

27.78 ± 1.42 

26.05 ± 0.24 

24.32 ± 0.08 

23.38 ±0.04 

22.78 ± 0.02 

21.63 ±0.00 

21.14 ±0.00 


3812 

3 

32 36.66 

-27 42 58.5 

>28.02 

26.05 ±0.15 

24.87 ± 0.09 

23.75 ± 0.04 

23.35 ± 0.02 

22.23 ±0.01 

21.74 ±0.00 


4148 

3 

32 44.97 

-27 43 09.1 

27.33 ± 0.48 

25.00 ±0.05 

24.43 ± 0.05 

24.09 ± 0.04 

23.82 ±0.03 

22.63 ±0.01 

21.83 ±0.00 


4173 

3 

32 41.24 

-27 43 09.7 

>27.57 

26.31 ±0.30 

25.24 ±0.19 

24.48 ±0.11 

23.91 ±0.07 

23.29 ±0.03 

22.65 ± 0.02 


4324 

3 

32 31.32 

-27 43 16.1 

>27.48 

26.23 ± 0.29 

24.65 ±0.12 

23.85 ± 0.07 

23.23 ± 0.04 

22.25 ±0.01 

21.77 ±0.01 


4327 

3 

31 59.71 

-27 43 15.5 

>28.40 

28.48 ± 1.29 

26.81 ±0.47 

26.89 ±0.61 

25.90 ±0.36 

25.47 ±0.17 

24.37 ± 0.08 


4534 

3 

32 29.99 

-27 43 22.7 

>28.07 

27.44 ±0.51 

25.42 ±0.13 

24.37 ± 0.06 

23.84 ± 0.04 

22.74 ±0.01 

22.26 ± 0.01 


4648 

3 

32 26.10 

-27 43 26.7 

28.39 ± 1.28 

26.60 ± 0.23 

25.47 ±0.13 

24.77 ± 0.08 

24.29 ± 0.07 

23.49 ± 0.02 

22.92 ± 0.02 


4846 

3 

32 25.93 

-27 43 31.1 

>28.28 

>28.37 

26.45 ±0.31 

26.23 ±0.31 

25.28 ±0.17 

24.08 ± 0.04 

23.67 ± 0.03 


4921 

3 

32 02.44 

-27 43 35.8 

>28.03 

>28.14 

>27.58 

26.98 ± 0.75 

27.88 ± 1.69 

26.13 ±0.23 

24.91 ±0.10 


5410 

3 

32 32.11 

-27 43 55.3 

>27.22 

26.33 ±0.38 

25.26 ±0.24 

25.83 ±0.50 

25.53 ±0.43 

23.87 ±0.06 

23.45 ± 0.05 


5529 

3 

32 01.77 

-27 44 01.1 

27.56 ±0.88 

28.02 ± 1.16 

25.76 ± 0.24 

24.76 ±0.12 

24.39 ± 0.09 

23.46 ± 0.03 

22.87 ± 0.02 


5685 

3 

32 12.78 

-27 44 07.7 

27.98 ± 1.78 

27.70 ± 1.10 

25.17 ±0.18 

24.95 ±0.18 

24.17 ±0.09 

23.44 ± 0.03 

23.06 ±0.03 


5735 

3 

32 02.83 

-27 44 09.7 

>28.05 

26.80 ± 0.29 

26.88 ±0.52 

26.52 ± 0.45 

26. 10 ±0.33 

25.80 ±0.17 

24.66 ± 0.08 


6845 

3 

32 06.57 

-27 45 14.0 

>27.57 

27.81 ± 1.15 

25.93 ± 0.34 

25.06 ±0.18 

24.80 ±0.15 

23.44 ± 0.03 

22.69 ± 0.02 


7202 

3 

32 06.40 

-27 45 54.7 

27.10 ±0.61 

25.76 ±0.15 

24.78 ±0.10 

23.67 ± 0.05 

23.16 ±0.03 

22.37 ±0.01 

21.71 ±0.01 

X-ray ID 82 b 


Notes. 

a Coordinates refer to the J2000 epoch. 
b X-ray identifications are from Luo et al. (2008). 


explicitly include the shot noise from the objects. We perform 
all fits with GalFit in units of counts using the H - band mosaic 
effective exposure time of 5017.7 s. 

We excise an 81 x 81 pixel 2 stamp centered on each galaxy; 
this size was chosen as a compromise between a sufficient 
number of sky pixels for a robust sky estimation by GalFit, 
and the number of (generally unassociated) nearby galaxies. 
Neighboring galaxies can bias the parameter estimation of the 
primary galaxy, therefore one must carefully mask out the 
unmodeled objects, or simultaneously fit all the objects in 
the stamp. We opt for the latter, since it avoids the ambiguities 
in pixel masking, and permits the flux in a given pixel to be 
represented as the sum of independent components. However, 
without good initial conditions, the GalFit algorithm may not 
converge to a meaningful solution, since the number of degrees 
of freedom can be very large. Therefore, to estimate better initial 
conditions for each galaxy in the stamp and ultimately ensure 
convergence in the final simultaneous solution, we fitted the 
two-dimensional light profiles of the primary BzH galaxy in 
question and any neighboring galaxies in a multi-stage process 
as follows. 

1. We identify all neighboring galaxies and their associated 
pixels with SExtractor, using the same settings men- 
tioned in Section 2.2. Any sources with FWHM_IMAGE# ^ 1 
are eliminated from the SExtractor segmentation maps 


and object catalogs, since they are likely cosmic rays or 
other artifacts. These objects are masked in all subsequent 
fitting processes. 

2. The pixels of any galaxy whose isophotes are truncated 
(based on the SExtractor flags) are masked, by setting 
the uncertainty maps to 10 10 ADU. These galaxies are no 
longer considered in the GalFit process. 

3. We model the light distribution of each remaining galaxy 
(including the primary BzH galaxy) individually, while 
masking all the pixels associated with every other galaxy. 
For any neighboring galaxy with semi-minor axis of 
ELIMAGE# ^ 1 pixel from SExtractor, we switch from 
fitting a Sersic profile to a PSF model, in order to eliminate 
degenerate degrees of freedom. Our results are robust to 
the choice of semi-minor axis limit, provided that we do 
not permit it to be larger than the size of the empirical PSF 
(discussed in more detail below). It is important to note that 
this step is only present to get reasonably accurate initial 
conditions for subsequent simultaneous object fitting. 

4. We refit a GalF it model, which contains a combination of 
point sources and Sersic profiles, to the stamp as a whole 
using the results from the previous step as the initial guesses. 


4 


In a few cases, we manually masked diffraction spikes or stellar 
halos clearly associated with foreground stars, which were just 
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Figure 2. BzH color images. Each stamp is « 4'.'5 38 kpc at z = 1 .6) on a side, has north up and east left, and has a pixel scale of 0?090 pixel -1 . All images are 

shown with the same color and logarithmic intensity scale. We find the object just to the south of 2750 (and 2749) is just barely too blue to have been included in 
our sample, which is consistent with the suggestion that this object has recently (~150 Myr) undergone an intense burst of star formation of 500-2000 Mq yr 
(van Dokkum & Brammer 2010; Ferreras et al. 2011). 


outside the field of view of the stamp before proceeding through 
the above procedure. 

For the individual and simultaneous fits, we placed con- 
straints on various GalF it parameters to prevent the algorithm 
from diverging into an unphysical regime. We constrained the 
centroid of any component to be within the ±2er of the cen- 
troid determined by SExtractor, and the total magnitude to be 
35 ^ H tot A- MAG AUTO/./ 2 mag. In general, these constraints are 
so weak that they generally play no role in the fitting whatsoever, 
but they reduce the sensitivity to the pixel masking with the seg- 
mentation maps. We additionally constrain the Sersic index and 
effective radius to be 0.01 ^ n ^ 8 and 0.01 ^ r e /A_IMAGE ^ 5, 
respectively. These constraints are considerably stronger, and 
we recognize that the model is likely incorrect when GalFit 
converges to a solution which is on these boundaries. These 
cases are rare, and generally a sign that additional astrophysi- 
cal components are needed (e.g., bulge/disk separately, nuclear 


point sources, or merger signatures), that the frame was inappro- 
priately sized, that the object was unresolved, and/or that there 
was some additional light component or image defect present in 
the image (e.g., diffraction spikes, stellar halos, and/or cosmic 
rays). 

As our primary interest here is on the sizes of these galaxies, 
it is imperative that we ensure the effective radii are robustly 
measured. While GalFit reliably determines the random un- 
certainty which follows from the maximum likelihood analysis, 
the total uncertainty should include a systematic term as well. 
To estimate the contribution from the systematic uncertainty, 
we construct a grid of simulated galaxies with brightnesses 
20 mag ^ H ^ 25 mag, effective radii 0.5 pixel sA r e <) 5.5 pixel, 
and a fixed Sersic index of n = 4 (for simulation only). These 
galaxies are convolved with the PSF and embedded in a blank 
region of the H - band mosaic. We then fitted these simulated 
galaxies with GalFit using the above procedure and find that 
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Figure 3. BzH galaxy counts. We show the differential surface density of 
galaxies selected based on the usual BzK color criteria (Daddi et al. 2004) and 
a typical color of (( H — K )) = 0.7 mag, assuming an instantaneous burst of 
star formation and passive evolution from zf orm =10. For comparison, we show 
the BzK object counts from Lane et al. (2007) from the UKIDSS/UDS survey 
as open symbols, which were derived from a strict BzK selection. The upper 
axis is simply a transformation using the assumed ((H — K)) color and the right 
axis applies to our objects. To demonstrate that the apparent plateau in the BzH 
counts // A n ~ 22 mag is not due to sample incompleteness in the H band, we 
show the total source counts from the ERS data as a dashed line (Windhorst 
et al. 2011). Moreover, this plateau is unlikely related to the z'-band data, since 
two of our selection criteria (H Sj 25 mag and (V — 11) 1.8 mag) restrict 

our sample to have z! F 26.8 mag, which is roughly 0.5 mag brighter than the 
z'-band limit (Dahlen et al. 2010). The uncertainties presented here reflect the 
Poisson variation in the object counts and do not include any contribution from 
cosmic variance. 


GalFit typically overestimates the effective radii by ~10% and 
underestimates the uncertainties by ~ 1 0%. Both the systematic 
bias and uncertainty are weakly brightness dependent, where 
fainter galaxies are more biased and uncertain. We remove the 
brightness-dependent bias and quadratically add the uncertainty 
to our effective radii measurements reported in Table 2. 

Many of the these red galaxies are very small and, even with 
the space-based imaging, may still be unresolved. Therefore, 
it is critical to properly identify which galaxies are resolved 
and have reliable effective radii measurements. We begin by 
swapping the fully variable Sersic model for the primary galaxy 
with a pure PSF model, which can only vary in position and 
brightness. However, as we have simultaneously fitted every 
object in our postage stamps, we anticipate that GalFit may 
incorrectly change the parameters of neighboring (unrelated) 
galaxies to compensate for the poor primary model, particularly 
in the case of a well-resolved primary galaxy. Therefore, when 
we use the PSF model for the primary galaxy, we hold the 
parameters of the neighboring sources fixed at the values found 
previously by GalFit. We now compare the goodness-of-fit 
statistics for these Sersic and PSF models by considering the 
quantity: 

P _ x PSF ^ Sersic ^ 

^ Searsic 

and expect that sources with low values of F are equally well 
characterized by a PSF model as by the more complex Sersic 
profile. We calibrate this quantity by computing the /’’-values of 
the known ERS stars, which were used to derive the empirical 
PSF used above. In Figure 4, we show the F - values as a function 
of the GalF it-derived PSF magnitude, with the galaxies plotted 
as filled-blue points and the stars as red asterisks. The trend is as 



Figure 4. Comparison of PSF and Sersic models. We show the fractional 
difference between the x 2 goodness-of-fit measures from GalFit for the PSF 
and Sersic models. The BzH galaxies are shown as filled-blue points, while 
the red asterisks represent 75 ERS point sources, presumed to be stars. For 
GalFit to optimally fit an unresolved source with a PSF, it will drive the 
effective radius and Sersic index to unphysical regimes, de facto fitting a PSF. 
Therefore, the goodness of fit of a Sersic model and a PSF should be roughly 
equal. However, if one were to fit a PSF to a resolved object, then there should 
be a noticeable increase in the goodness-of-fit statistic. We adopt a critical value 
of F cr j t = 0.025, see Section 4.1 for more details. Of the 75 stars used in this 
test, only 4 have F > F cr i t suggesting this is a reliable way to classify unresolved 
objects. 

(A color version of this figure is available in the online journal.) 


expected: nearly all point sources can be found at — 0.5 < F < 0, 
while the BzH galaxies are generally at F > 0 which depends on 
brightness. Therefore, we define objects which can be equally 
characterized by a PSF as by a Sersic fit as having F F :: F crit , 
while objects with F > F a y are more extended than the known 
ERS stars. We adopt F cnl — 0.025 and note that only 4/75 stars 
have F > F cnu which is consistent with Bond et al. (2009). We 
give the GalF it results in Table 2. 

4.2. Stellar Populations and Photometric Redshifts 

We fit the matched-aperture HST photometry with a library 
of stellar population synthesis models to simultaneously deter- 
mine the stellar mass, population age, and redshift with our own 
software. Our model grid consists of a four-dimensional pa- 
rameter space spanned by redshift (z), stellar population age 24 
(/), V-band extinction (Ay), and star formation timescale (r) for 
an exponentially declining star formation history. We assume 
solar metallicity, a Salpeter initial mass function (IMF), and 
adopt the CB07 population synthesis models. The allowed pho- 
tometric redshifts ranged from 0 ^ z < 7 with A z = 0.01, the 
ages presented by CB07, extinctions of Omag ^ Ay ^ 2 mag 
with A A v = 0.2 mag, and star formation timescales of r e 
[10 -3 , 1, 2, 3, 5, 15, 30, 10 3 ] Gyr. We include a 10% systematic 
uncertainty on the observed fluxes to account for uncertainties 
in the zero points and assumed templates. We compute the la 
uncertainties on the stellar population parameters (e.g., mass, 
age, etc.) by use of a simple Monte Carlo calculation. For each 
band for a given galaxy, we draw a normal random variable 
with mean and standard deviation equal to the flux and flux 
uncertainty (without the systematic term) and catalog the new 
population parameters. By repeating for 10 4 iterations, we build 


24 We impose the usual self-consistency constraint that the a galaxy be 
younger than the age of the universe. 
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ID 

r a 
' e 

(arcsec) 

n b 

Xv 

F° 

Notes 

408 

0.04 ±0.01 

4.72 ± 0.45 

0.487 

0.160 


606 

0.09 ±0.01 

6.92 ± 0.92 

0.404 

0.158 


1696 

0.44 ± 0.48 

1.21 ±0.10 

0.356 

0.170 


2227 

0.36 ± 0.04 

2.93 ± 0.09 

0.364 

2.517 


2377 

0.18 ±0.05 

4.17 ± 1.01 

0.322 

0.049 


2749 

0.51 ±0.19 

0.57 ± 0.05 

1.849 

0.191 


2750 

1.02 ±0.26 

5.65 ±0.58 

1.891 

0.591 


2871 

0.03 ± 0.01 

8.00 ± 2.76 

0.567 

0.023 

Unresolved 

3000 

0.60 ±0.06 

8.00 ±0.14 

0.814 

9.440 

Tidal tail — the fit is dubious 

3152 

0.14 ±0.01 

4.66 ±0.12 

0.515 

3.323 


3237 

0.22 ± 0.02 

7.54 ±0.30 

0.432 

2.268 


3360 

0.11 ±0.01 

7.58 ± 0.46 

0.466 

1.382 


3376 

0.17 ±0.05 

2.21 ±0.44 

0.349 

0.054 


3471 

0.08 ±0.01 

4.40 ±0.16 

0.478 

1.525 


3488 

2.85 ± 6.07 

0.73 ±0.13 

0.494 

0.058 


3551 

0.25 ± 0.03 

8.00 ±0.22 

0.619 

5.687 

Additional nuclear point source was included 

3812 

0.16 ±0.02 

3.12 ±0.07 

0.486 

3.767 


4148 

0.04 ± 0.01 

5.73 ±0.32 

0.443 

0.690 


4173 

0.38 ±0.05 

0.98 ± 0.02 

0.359 

2.373 


4324 

0.25 ± 0.02 

3.27 ± 0.08 

0.404 

4.536 


4327 

0.03 ± 0.02 

2.96 ± 3.62 

0.411 

0.000 

Unresolved 

4534 

0.14 ±0.01 

2.21 ± 0.07 

0.388 

2.718 


4648 

0.18 ±0.02 

2.08 ±0.10 

0.355 

0.913 


4846 

0.17 ±0.02 

1.30 ±0.13 

0.506 

0.232 


4921 

0.18 ± 0.10 

8.00 ±4.50 

0.337 

0.012 

Unresolved 

5410 


2.64 ± 1.00 

0.364 

-0.009 

Likely three distinct clumps — the fit is dubious 

5529 

0.28 ± 0.03 

2.63 ±0.12 

0.369 

1.107 


5685 

0.47 ±0.21 

0.36 ±0.02 

0.386 

1.399 


5735 

0.09 ± 0.02 

8.00 ±4.15 

0.349 

0.012 


6845 

0.43 ± 0.09 

1.33 ±0.04 

0.385 

1.648 


7202 

0.30 ±0.03 

2.29 ± 0.04 

0.434 

5.825 



Notes. 

a The uncertainties on the effective radius include the systematic uncertainties determined in Section 4. 1 . 
b The Sersic index in Equation (3). 

c The fractional difference between the goodness of fit for the PSF and Sersic models (see Section 4.1 for more details). 


up a distribution of each stellar population parameter and take 
the mean and standard deviation as the measured quantity and 
uncertainties, respectively. 

We apply this approach to the complete sample of 30 
BzH galaxies. In Figure 5, we show the photometric redshift 
distribution and the comparison to a spectroscopic redshift, 
when available. As expected, the median redshift of the sample 
is Zphot = 1 .6 ±0.6, 25 where the uncertainty reflects the standard 
deviation of the distribution. To estimate our typical redshift 
uncertainty, we compare to the published spectroscopic redshifts 
in the lower panel of Figure 5. Based on the root-mean- squared 
scatter of the seven objects with known spectroscopic redshifts, 
we estimate our uncertainties in ( 1 -Hz) are ~3.3%. We do recover 
many of the known galaxy cluster members at z ~ 1.6 (Kurk 
et al. 2009). We present our resulting photometric redshifts and 
stellar population parameters in Table 3. 

5. BzH GALAXY SIZE EVOLUTION 

Based on the photometric redshift estimates and the spec- 
troscopic data (where available), the BzH selection reliably 
identifies galaxies in the interval z p hot~ L6. However, to study 



I , . . . I , . JjL I ± , , , I , , , , I , , , , 

0.5 1.0 1.5 2.0 2.5 3.0 

Photometric Redshift 


Figure 5. Photometric redshift distribution. In the upper panel, we show 
the photometric redshift distribution derived by the procedure described in 
Section 4.2. As expected, our sample of BzH galaxies is generally located at 
Zphot = 1.6. In the lower panel, we show the fractional difference between our 
redshift estimates and spectroscopic redshifts (where available) as a function of 
their photometric redshift. Based on the rms on the fractional differences, we 
estimate that (1 + z) is accurate to ~3.3%. 


25 This could be potentially biased given the presence of a known overdensity 
at 1.6 (Kurk etal. 2009). 
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Table 3 

Photometric Redshifts and Stellar Population Parameters 


ID 

£phot 

£spec 

log Age a 
(log yr) 

log Mass b 
(log Mq) 

Xv 

408 

1.58 ±0.05 


8.7 ±0.2 

10.6 ±0.1 

0.9 

606 

1.89 ±0.19 


8.8 ±0.3 

10.6 ±0.2 

1.9 

1696 

3.78 ± 0.83 


8.7 ±0.2 

11.4 ±0.7 

3.3 

2227 

0.79 ± 0.44 


9.1 ± 0.0 

9.9 ±0.6 

4.8 

2377 

2.35 ± 0.20 


8.7 ±0.3 

10.5 ± 0.7 

3.0 

2749 

1.96 ±0.20 


8.4 ± 0.5 

10.9 ± 0.2 

4.2 

2750 

1.90 ±0.10 


8.2 ±0.3 

11.1 ±0.2 

6.8 

2871 

1 .52 ± 0.05 


8.7 ±0.2 

10.4 ±0.1 

2.2 

3000 

1 .50 ± 0.05 


8.9 ±0.1 

11.5 ± 0.1 

0.6 

3152 

1.51 ±0.06 

1.367 

9.0 ± 0.1 

11.2 ±0.1 

1.4 

3237 

1 .63 ± 0.07 

1.615 

8.8 ±0.3 

11.1 ± 0.1 

1.9 

3360 

1 .73 ± 0.20 


8.9 ±0.1 

11.1 ± 0.1 

0.9 

3376 

1 .48 ± 0.09 


8.5 ± 0.3 

9.8 ± 0.2 

0.9 

3471 

1 .64 ± 0.05 

1.610 

8.9 ±0.1 

10.9 ±0.1 

1.0 

3488 

4.73 ± 0.29 


8.7 ±0.1 

11.4 ±0.4 

2.5 

3551 

1.58 ±0.06 


9.0 ± 0.1 

11.3 ± 0.1 

1.7 

3812 

1.71 ±0.12 

1.614 

9.0 ± 0.1 

11.0 ±0.1 

1.7 

4148 

2.36 ± 0.04 


8.6 ±0.1 

11.0 ±0.1 

2.9 

4173 

0.98 ± 0.42 


9.0 ± 0.3 

10.1 ± 0.7 

1.2 

4324 

1.54 ±0.08 


9.0 ± 0.1 

10.9 ±0.1 

2.7 

4327 

3.13 ±0.76 


8.6 ±0.3 

10.9 ±0.8 

1.8 

4534 

1.53 ±0.10 

1.604 

9.0 ± 0.1 

10.8 ± 0.1 

1.8 

4648 

0.86 ± 0.62 


8.8 ±0.4 

9.7 ±0.8 

2.5 

4846 

1.63 ±0.10 


8.8 ±0.2 

10.4 ±0.2 

3.5 

4921 

3.57 ± 0.88 


8.8 ±0.2 

11.4 ± 1.0 

2.6 

5410 

1.99 ±0.18 


8.5 ± 0.3 

10.4 ± 0.3 

4.7 

5529 

1.17 ±0.41 


9.0 ± 0.1 

10.4 ± 0.2 

1.0 

5685 

1 .48 ± 0.07 


9.0 ± 0.1 

10.3 ± 0.2 

0.9 

5735 

2.66 ±0.12 


8.4 ±0.1 

9.9 ±0.2 

1.0 

6845 

2.03 ± 0.53 


8.9 ±0.2 

10.8 ± 0.2 

1.0 

7202 

1.45 ±0.11 

1.329 

8.4 ±0.1 

11.0 ±0.2 

1.2 


Notes. 

a Population age assuming an exponential star formation history. 
b Stellar mass. 


the evolution of their sizes with redshift, we must compare to 
other similarly selected samples. The high-redshift (z > 1.5) 
samples are generally derived from similar color criteria pre- 
sented here and have effective radii measured in the H band, for 
a rest-frame wavelength of A rest ~ 6500 A. To ensure fair com- 
parisons with lower redshift samples, we require sizes measured 
at similar rest-frame wavelengths. For the low-redshift data, we 
use the sample of 8666 early-type galaxies at z ~ 0.2 with 
effective radii measured in the i' band (Bernardi et al. 2003) 
selected from the Sloan Digital Sky Survey, Early Data Re- 
lease (SDSS-EDR; York et al. 2000). The stellar masses for the 
7th Data Release 2 ' 1 galaxies were determined following Salim 
et al. (2007). By cross matching these samples, we obtain 8595 
galaxies suitable for our low-redshift comparison. We select 
three mass ranges which are volume limited based on the upper 
and lower flux limits imposed by detector saturation in SDSS 
and the Bernardi et al. (2003) brightness criterion, respectively. 
In Figure 6, we show these volume-limited selections (black 
boxes) and the flux limits (dashed lines). If these limits are not 
strictly imposed, then an artificial redshift dependence on the ef- 
fective radii will be introduced as the radii are tightly correlated 
with the stellar masses, which roughly scale with luminosity. 
For example, without these volume limits, the effective radii of 
the SDSS galaxies will seem to increase with redshift since the 

26 Obtained from http://www.mpa-garching.mpg.de/SDSS/DR7/. 
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Figure 6. Selection of SDSS galaxy sample. We show the stellar mass as a 
function of spectroscopic redshift for the early-type galaxy sample of Bernardi 
et al. (2003) drawn from the SDSS-EDR. The stellar mass estimates were derived 
from stellar population fits to the SDSS photometry (Salim et al. 2007) and are 
similar in nature to those described in Section 4.2. The dashed lines indicate 
the bright (i f 14 mag) and faint ( i' ^ 17 mag) completeness limits set by 
detector saturation and the Bernardi et al. (2003) selection, respectively. The 
boxes show our volume-limited selections for proper low-redshift comparisons. 
If such limits are not imposed, then an artificial trend in effective radius with 
redshift will arise from a Malmquist-type bias. 



survey is not sensitive to the lower mass (smaller) galaxies at the 
higher redshifts. Eliminating this critical Malmquist-like bias in 
the low-redshift sample ensures a fair comparison between the 
high- and low-redshift data. 

To investigate the passively evolving galaxy size evolution 
at a fixed stellar mass, we show in Figure 7 the effective radii 
as a function of redshift for this work (filled circles) and a 
host of published measurements as described in the caption. We 
convert all mass estimates from the literature to the Salpeter 
IMF based on corrections discussed in the Appendix. For clar- 
ity, we do not show three objects whose measured effective 
radii are r e ^ 0.5 pixel as marked by the gray shaded re- 
gion. While these objects are formally consistent with being 
extended, such small sizes are dubious and highly degenerate 
with the Sersic index and local sky estimation. Furthermore, 
there are 1 1 objects from our sample which are either outside 
our mass range (10.3 log (M*/M 0 ) ^ 11.5) or have photo- 
metric redshifts at z p h 0 t ^ 2.6, as the spectral energy distribution 
fits are highly suspect. We overplot the power-law models of the 
form R e oc (1 + z)~ a (thick black line) where a is determined 
for each mass range separately and present the goodness-of- 
fit statistic for the two models in the lower-left corner of each 
panel. We treat the canon of SDSS as a single low-redshift 
constraint, using their mean and error in the mean for the fits. 
For comparison, we show the mass-independent model (dashed 
blue line) in each panel, which was determined by fitting the 
data in all three panels simultaneously. We note that fits of the 
form R e oc H (z) _/i (where H(z) is the Hubble parameter) 
give a similar qualitative result: the amount by which galaxies 
are smaller in the past depends on their stellar mass. Newman 
et al. (2010) identified a similar result for 12 galaxies with dy- 
namical masses determined from measured velocity dispersions. 

Our preliminary estimate of a mass-dependent size evolution 
model is qualitatively similar to the predictions from the semi- 
analytic models of Hopkins et al. (2009), where the power- 
law index varies with mass as a(M*)« 0.23 log (M*/10 9 M Q ). 
However, their model underpredicts our measured power-law 
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Figure 7. BzH galaxy size evolution. We show the effective radius vs. redshift 
for the various stellar mass slices described in Figure 6 and several early- 
type galaxy samples: SDSS sample (Bernardi et al. 2003; Salim et al. 2007; 
small dots), Daddi et al. (2005; asterisks), Longhetti et al. (2007; diamonds), 
Toft et al. (2007; downward triangles), Zirm et al. (2007; closed squares), 
Cimatti et al. (2008; open squares), Damjanov et al. (2009; upward triangles), 
Carrasco et al. (2010; open circles), and this work as large-filled circles. All 
stellar masses assume a Salpeter IMF with solar metallicity (see Section 5 
and the Appendix for the conversion factors). We show the power-law fit 
R e oc (1 + z)~ a , where a is determined separately for each mass bin (thick 
black line) and for all bins together (dashed blue line), and report the reduced 
X 2 for each model in the lower left of each panel. The gray shaded region in 
each panel indicates the position of unresolved objects with r e <0.5 pixel; 
we omit three objects which are marginally resolved (based on their F-ratios), 
but are within this region. The red stars indicate the average and rms values in 
three redshift intervals, which are only used to illustrate the trend. Based on 
the reduced x 2 values, we conclude the mass-dependent model (full black line) 
better reproduces the relationship between the effective radius and redshift than 
the mass-independent (dashed blue line), while preserving the local R e -M * 
trend (e.g., Shen et al. 2003). Since the majority of high-redshift objects have 
very similar effective radii, stellar mass, and redshift distributions, much of the 
variance in a is set by the SDSS sample. However, larger high-redshift samples 
(comparable in size to those at low redshift) are needed to assess their relative 
importance in setting the mass dependence on a. 

(A color version of this figure is available in the online journal.) 

indices in our three mass regimes, but does give the same 
qualitative steepening of a with mass. Following their approach, 
we derive tentative relationships for the power-law indices: 

odM,) * (-0.6 ± 0.7) + (0.9 ± 0.4) x log ( (5) 

and caution that with only three independent mass bins, this 
scaling should be considered preliminary at best. However, 
this result gives indices consistent with the reported value of 
Buitrago et al. (2008), who find that a(M* > 10 11 M e ) & 
1.5. In Section 7, we discuss some potential shortcomings 
of the current data and propose experiments for future data 
sets. 

6. DISCUSSION 

It has become relatively well established that early-type 
galaxies are indeed smaller at high redshift than their local 


counterparts for a given stellar or dynamical mass (e.g., 
Trujillo et al. 2009; Newman et al. 2010). However, the causes 
for this result are far less clear, or agreed upon. The likely mech- 
anisms can be broadly characterized as astrophysical effects 
(early-time major mergers, late-time minor mergers, adiabatic 
expansion, or stellar mass-to-light gradients), or observational 
biases (underestimating the effective radii, overestimating the 
stellar masses, or incorrectly assuming that the high- and low- 
redshift populations are directly comparable). Using a semi- 
analytical model, Hopkins et al. (2010) find that the factor of 
~5 increase in the effective radii over the last ~10 Gyr for 
galaxies with M ; , 1 0 1 M Q can be explained by a combina- 

tion of these effects, with the late-time minor merging play- 
ing the largest role. However, we find that lower mass sys- 
tems exhibit notably weaker redshift evolution, suggesting a 
different mixture of the main processes may be at work. For 
example, the adiabatic expansion mode may become more crit- 
ical given the shallower gravitational potentials in the low-mass 
systems. 

As noted above, the tentative stellar inass-dependent size 
evolution presented in this work is similar to the findings of 
Newman et al. (2010) for dynamical masses estimated from 
velocity dispersions, and has an interesting consequence for the 
R e -M * relation. Shen et al. (2003) find that SDSS elliptical 
galaxies follow the scaling relation R r , = R e \\{M*/\0 n M 0 ) y , 
where R^n = 4.16 kpc and y = 0.56. However, the mass- 
dependent size evolution found by us and by Newman et al. 
(2010) imply either a fundamental change in the R e -M scaling 
relationship at high redshift (such as a flattening of the effective 
radius for decreasing stellar mass) or a redshift-dependent value 
of y , in addition to the usual lower value of R e ,n- Unfortunately, 
at this stage the high-redshift data cannot distinguish these two 
scenarios or shed light on the cause for this change. In fact, given 
the large scatter in R e -M * at high redshift (z~ 1.5), most of the 
size evolution is driven by the Bernardi et al. (2003) sample at 
z-0.2. 

The favored red galaxy formation paradigm suggests that 
passively evolving galaxies have come from a population of 
more active systems in their recent past. These systems originate 
from mergers of gas-rich disks that trigger intense starbursts 
which are later quenched by active galactic nuclei (AGNs), 
finally resulting in dead spheroidal systems spent of their gas and 
no longer forming stars (e.g., Faber et al. 2007, and references 
therein). The Lyman break galaxies (LBGs) observed at z > 3 
are a possible progenitor system (e.g., Overzier et al. 2010). 
While their stellar masses are far less well constrained due 
to the lack of rest-frame IR data, the majority of LBGs have 
M* ~ 10 10 M e between 3 < z < 6 (e.g., Papovich et al. 2001; 
Yan et al. 2006). Therefore, the typical LBG at 3 < z < 6 
belongs to the lower panel of Figure 7 and has effective radii 
of R e ~ 1 kpc (e.g., Ferguson et al. 2004; Bouwens et al. 
2006; Hathi et al. 2008; Oesch et al. 2010), after transforming 
to an equivalent rest-frame wavelength (Barden et al. 2005). 
The LBGs are then consistent with our size evolution model, 
possibly suggesting that similar physical mechanisms driving 
the passive galaxy evolution may be at work for the LBGs 
as well. Given their increased star formation rates, larger gas 
content, and lower total mass, the mass-loss modes, whether 
driven by AGNs (e.g.. Fan et al. 2008) or stellar winds (e.g., 
Damjanov et al. 2009), are potentially more important. However, 
we recognize that these samples (LBGs and passive galaxies) 
are quite different and that the next generation of space-based 
infrared instruments (Near-Infrared Camera and Mid-Infrared 
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Table 4 

Stellar Mass Transformations 3 


Ryan et al. 









log M*,i 






35c 

3 5 s 

3 6c 

3 6s 

37c 

37s 

75c 

75s 

76c 

7 6s 

77c 

log Mtj 

35s 

-0.25 












36c 

+0.06 

+0.31 











3 6s 

-0.19 

+0.06 

-0.24 










37c 

+0.11 

+0.36 

+0.05 

+0.29 









37s 

-0.14 

+0.11 

-0.20 

+0.04 

-0.24 








75c 

+0.07 

+0.32 

+0.01 

+0.25 

-0.04 

+0.20 







75s 

-0.17 

+0.08 

-0.24 

+0.01 

-0.29 

-0.04 

-0.25 






76c 

+0.06 

+0.31 

-0.00 

+0.24 

-0.06 

+0.19 

-0.01 

+0.24 





7 6s 

-0.18 

+0.07 

-0.25 

-0.00 

-0.30 

-0.05 

-0.25 

-0.01 

-0.24 




lie 

+0.09 

+0.34 

+0.02 

+0.27 

-0.01 

+0.24 

+0.02 

+0.27 

+0.03 

+0.27 



77s 

-0.16 

+0.09 

-0.22 

+0.03 

-0.25 

-0.01 

-0.22 

+0.02 

-0.22 

+0.03 

-0.24 


Note. 

a We report the quantity | log 


Instrument on James Webb Space Telescope , JWST) will provide 
a much clearer picture for the high-redshift (z > 3) size 
comparisons. 

7. SUMMARY 

We identified 30 passively evolving galaxies to H ^ 25 mag 
from a set of color criteria similar to the p BzK selection. We 
measure rest-frame optical (/. rest ~ 6500 A) effective radii as 
two-dimensional fits to the //-band image. By comparing with 
several other comparable samples at similar redshifts, we find 
that the size evolution depends on the stellar mass. We give 
tentative scalings between the power-law index of the R e -z 
relation and stellar mass. While we have made every attempt 
at matching the selection biases between the high- and low- 
redshift samples, subtle differences may remain (such as the 
requirements on spectroscopic redshift and visual classification 
on the low-redshift data not matched at high redshift). Given the 
various requirements on placing SDSS fibers and the seeing- 
limited imaging, it is easy to believe that there may be systematic 
differences in the way the samples are created. Using future 
surveys with ultraviolet imaging with HST (such as planned 
observations with CANDELS and of the Hubble Ultra-Deep 
Field), we can extend the redshift range to 0.3 < z < 1.3 while 
consistently selecting galaxies of identical rest-frame spectral 
properties and removing any potential differences between the 
high- and low-redshift samples. 
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APPENDIX 

STELLAR MASS CONVERSIONS 

In most cases when authors determine stellar masses from 
population synthesis techniques, they adopt a single IMF and 
metallicity with the understanding that the stellar masses are 
dependent on these choices. Since there is no clear indication 
of which IMF may be appropriate, and modeling an arbitrary 
IMF and metallicity with broadband observations is prone to 
large degeneracies (Conroy et al. 2009), we are left to convert 
all stellar masses taken from the literature to a common set of 
assumptions. To derive these conversions, we refit our 30 objects 
with a different IMF and metallicity and compute 

= (A1) 

where M* ,■ and M* j are the resulting stellar masses for the 
different input models i and j , respectively. We extend this 
analysis to include the Bruzual & Chariot (2003, BC03) models, 
which did not include flux contributions for thermally pulsating 
asymptotic giant branch stars. In Table 4, we tabulate the C 
matrix to convert the mass estimates obtained from the literature 
(as discussed in Section 5), where a model is encoded by 
three characters: the first indicates the vintage (3 = BC03 and 
7 = CB07), the second indicates the metallicity (5 = 0.4 Z Q , 
6 = 1.0 Z Q , and 7 = 2.5 Z Q ), and the third indicates the IMF 
(S = Salpeter or C = Chabrier). 
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